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Abstract 

Chiral symmetry is known to be decisive for an understanding of the 
low energy sector of strong interactions. It is thus important for a 
model of relativistic heavy ion collisions to incorporate the dynamical 
breaking and restoration of chiral symmetry. Thus we study an expan- 
sion scenario for a quark-meson plasma using the Nambu-Jona-Lasinio 
(NJL) model in its three flavor version. The equations of motion for 
light and strange quarks as well as for pions, kaons and etas are solved 
using a QMD type algorithm, which is based on a parametrization of 
the Wigner function. The scattering processes incorporated into this 
calculation are of the types qq <-> qq, qq <-> qq, qq <-> MM and M <-> qq. 



*To appear in the proceedings of the 4th International Symposium on Strangeness in 
Quark Matter, Padova, Italy, July 20-24. 
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In relativistic heavy ion collisions, a very hot and dense system of strongly in- 
teracting matter is created. At these high temperatures and densities, the state of 
matter changes drastically. The two main effects which are predicted from quantum 
chromodynamics (QCD) are (i) deconffnement and (ii) chiral symmetry restaura- 
tion. Lattice simulations indicate, that these two effects are not independent of 
each other, but occur at the same time ]l| . Although experiments searching for this 
new state of matter have been performed at the Cern SPS and will be performed at 
RHIC and LHC, the theoretical description still leaves many open questions. Some 
of these are (i) the question, if a local thermal and/or chemical equilibrium will be 
established, (ii) the influence of the in-medium changes of hadronic properties and 
(iii) the description of the hadronization. The treatment of these problems is by 
no means simple, since a satisfactory, simple model for the description of the phase 
transition to a quark-gluon plasma (QGP) does not yet exist. Since the QGP, if 
created, is a transient state, a realistic model for its evolution should be able to 
handle nonequilibrium effects. A nonequilibrium theory of nonabelian gauge theo- 
ries such as QCD, however, cannot be provided using present days knowledge, since 
these theories even in thermal equilibrium are not yet sufficiently understood. Lat- 
tice gauge theories, on the other hand, have been successfully applied in order to 
extract the thermal behaviour of strongly interacting matter. They work, however, 
only in equilibrium situations. Another practical approach which is frequently used 
is hydrodynamics. These models, however, suffer from the fact that they work in 
the limit of infinitely many collisions per time interval, which might not be true in 
the first stages and surely is not true in the late stages of the collision. 

The ansatz presented here thus the modelization of the expansion of a hot system 
using an effective Lagrangian. The model interaction we use is the three flavor 
Nambu-Jona-Lasinio (NJL) model ||, defined by 



In this model, the exchange interaction has been contracted to give a pointlike 
interaction in flavor space. This approximation works as long as the effective mass 
of the gluon is high. Recent QCD calculations support this assumption even at 
high temperatures ||. The most important feature of the Lagrangian (|]) is that 
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it preserves the chiral symmetry of QCD, i.e. in the limit m / — ^ it is invariant 
under transformations of the form 

ip -> exp (-iOaKlb) 4> ■ (2) 

It is a well established fact, that this symmetry is the essential ingredient for the 
description of the low lying hadronic states which are produced copiously in a heavy 
ion collision jp. 

The equilibrium properties of the NJL model in thermal equilibrium have been 
studied in great detail elsewhere so that only those properties will be men- 

tioned here, which are important for the understanding of the present article. For 
simplicity, first the case m / = will be detained. At temperature T = 0, chiral 
symmetry is spontaneously broken, which leads to a finite effective quark mass. As 
a consequence of the Goldstone theorem, Nj — 1 massless modes appear as quark- 
antiquark bound states. These states are identified with the 7r, K and rj mesons. 
At a certain finite temperature, however, chiral symmetry gets restored and one has 
again massless quarks, whereas the mesons become massive. In this phase, they 
do no longer exist as bound states but rather as resonances. If one introduces a 
finite current quark mass, m / 0, this picture gets slightly changed in that chiral 
symmetry is no longer an exact, but rather an approximate symmetry. Thus mesons 
have a small but finite mass at T = 0, as they have in nature, whereas the quark 
mass does not go to zero at large temperatures, but stays finite and for T — > oo goes 
to the current quark mass. 

The quark mass spectrum for this case is shown in Fig. The light quarks u 
and d are taken to be degenerate and will be denoted by the generic index q in 
the following. The current quark masses used in Fig. p] are mo q = 5.5 MeV and 
m 0s = 140 MeV. Due to chiral symmetry breaking, one has at T = an effective 
mass of m q = 368 MeV and m s = 550 MeV. In the temperature region T » 200 MeV, 
these masses drop and the effective mass of the light quarks becomes low. For 
the strange quarks, however, the situation is different. Due to the high current 
quark mass, the effective mass stays relatively high. At T = 350 MeV, which is 
far beyond any temperature to be realistically expected in present days heavy ion 
experiments, one still has m s = 300 MeV, which is about twice as high as the current 
quark mass. The pseudoscalar mesonic mass spectrum using the same parameters is 
shown in Fig. 0. At T = 0, the mesons have their physical masses of m n = 135 MeV, 
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mjf = 498 MeV and = 515 MeV. With rising temperatures, these values stay 
more or less constant, until to a certain temperature, where their masses become 
equal to the masses of their constituents, i.e. m n = 2m q , itik = fn q + m s or 
m n = 2m q , respectively. This happens when the m n and lines of Fig. |2| cross the 
2m q line, or the line crosses the m q + m s line, respectively, as is indicated by the 
arrows. At these temperatures, the mesons become unstable via a Mott transition 
0. At higher temperatures, they are unstable with respect to the decay M —>■ qq 
and obtain a finite width. The qualitative form of the mass spectrum of Fig. |2| has 
been confirmed by lattice calculations M. 

A finite temperature study of this model will surely not lead to an event genera- 
tor, since the model still contains free quark states, which are not observed in nature. 
It will, however, be able to study qualitative features of the plasma expansion, such 
as expansion time scales, the approach to equilibrium and thus the applicability of 
hydrodynamics, the production mechanism of hadrons etc. While the simulation 
program presently is limited to zero baryochemical potential, future extensions will 
remove this limitation and also allow for the study of strangeness destination and 
DCC formation. 

The treatment of the NJL model in nonequilibrium starts from the observation, 
that both quark and meson degrees of freedom can be simultaneously described by 
an equation of the Boltzmann type ||, 

(d t + d p Ed x - d x Ed p ) n{t, x, p) = J coU [n(t, x, p)} , (3) 

where ic n[ra] is a collision integral. In principle, one is thus able to describe a 
transition from a pure quark regime to a hadronic regime, where quarks are converted 
to hadrons via collision processes. The solution method we choose for Eq. is 



an algorithm of the QMD type |10[ , i. e. we parametrize the particle distribution 
functions by 



n(t,x,p) = Eexp ( - (f 2 *J t))2 ) exp (-^(p -p^f ) ■ (4) 



The center points of the distribution move on the characteristics of Eq. (|J), i.e. 

— * 

Xi(t) = Pi{t) I 'E Piif) — ~d x E + collision contributions . (5) 

Equation (151) has to be solved together with the gap equation 
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where the indices i, j and k run over all three quark flavors. For the computation 
of mesonic properties we take a shortcut by defining an effective temperature via 



using the equilibrium expressions, which are functions of T e fj, m q and m s . The 
collision processes entering Eq. (|5]) are (i) quark quark scattering qq <-> qq, qq <-> qq 
and gg <-> gg J7j, (ii) hadronization gg <-> MM ||, and (iii) meson decay M —>■ qq 
||. The latter process is only possible in the early phases of the expansion, when 
the effective temperature is higher than the Mott temperature. 

The initial conditions chosen presently do not contain strange quarks, whereas 
light quarks are distributed thermally within a sphere of a given radius. This kind 
of initial conditions has the immediate consequence, that the strange quark mass 
is even higher than in thermal equilibrium. The reason for this can be seen by 
writing Eq. (^) explicitly for strange quarks. One recognizes that the second term 
on the right hand side, which is proportional to G, only couples to the strange 
quark condensate. Since initially no strange quarks are present, this term does not 
receive medium corrections and thus makes the effective quark mass higher than in 
thermal equilibrium. Medium corrections do arise from the third term of Eq. (|J), 
which is proportional to the product of the up and down quark condensate. This 
contribution is, however, not sufficient to produce a large mass drop. 

Quantitatively, this can be seen from Fig. [| which shows the quark masses and 
the effective temperatures as a function of r at various times during the initial phase 
of the expansion. At t = 0, the light quark mass in the center of the system is low, 
according to the high particle density here. At larger radii, the mass goes up due to 
the gaussian form of the particle distribution. Note that the quark masses are only 
known at those points where a particle is present, thus the curves stop at the edge 
of the fireball. The effective temperature in the centre amounts to approximately 
250 MeV. At this temperature in equilibrium, one would expect a strange quark 
mass of 380 MeV according to Fig. |l|. In reality, however, one has a strange quark 
mass of approximately 450 MeV in the centre due to nonequilibrium effects. At later 



m q (x,t) = m q q (T e g(x,t)), where m q q (T) is the equilibrium form of the temperature 
dependence of m q as shown in Fig. 0. The mesonic properties can then be obtained 



5 



times, the system expands and thus the quark mass curve becomes flatter, while its 
width grows. Accordingly, the effective temperature drops. In the final state, the 
quark masses will tend towards the vacuum value, which means that the mean field 



part of the interaction ceases |TI 



Mesons are created during the evolution by collisions of quarks and antiquarks. 
Most of all mesons created are pions, which are most frequently produced by the 
collision of two light quarks. Kaons, on the other hand, are produced most easily 
by the collision of a light with a strange quark or two strange quarks, while the 
hadronization cross section for the creation of a kaon pair from a light quark an- 
tiquark pair is rather low ||. Eta mesons are most frequently created collisions of 
two light quarks, forming a pion and an eta. The time dependence of the multiplic- 
ities is shown in Fig. |j. It can be seen here, that the production of mesons starts 
immediately after the beginning of the expansion. The production rate is maximal 
at t = 0, when the particle density is high, giving thus rise to a large number of 
collisions. At later times the density drops and thus the number of collisions per 
time decreases. This in turn leads to a flattening of the multiplicity curves. 

Figure |5] shows the particle density, averaged over all solid angles, at time t = 
lOfm/c. It is clearly visible that the meson density is maximal at the same places 
where also the quark density is high. This confirms the picture gained from Fig. ^, 
i. e. that mesons are created within the bulk of the fireball. More insight can be 
gained from Fig. g, which shows the distribution of the effective temperatures at 
the creation points of the mesons for pions, kaons and etas respectively. All three 
distributions agree more or less with each other, up to the effects of lower statistics 
for the kaon and the eta. It can be seen that mesons are created most likely at 
temperatures directly below the Mott temperature, and within a temperature range 
of 50 MeV. This plot should be compared to the mean hadronization times for light 
and strange quarks in equilibrium, as have been calculated in M. It has been shown 
there, that these hadronization times have a minimum in the temperature range 
150MeV< T <200MeV, which agrees nicely with Fig. g. 

Figure [7] shows the density along the coordinate axes for t = and t = 30fm/c. 
It can be seen here, that, although one starts with a system, which is to a good ap- 
proximation spherically symmetric, one ends up with a final state which shows large 
fluctuations of the density with respect to the direction. This behaviour contradicts 
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hydrodynamics, where an initially symmetric system stays symmetric. 

To conclude, it hat been demonstrated how a chirally symmetric quark- meson 
plasma containing strangeness behaves out of equilibrium and how mesons are pro- 
duced. This investigation has been performed at zero baryochemical potential. Fu- 
ture investigations will remove this limitation and thus be able to study the mech- 
anism of strangeness destination. Also the study of DCC formation is planned. 
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FIGURES 

FIG. 1. Quark masses at finite temperature in equilibrium. Solid line: light quarks, 
dashed line: strange quarks. 

FIG. 2. Meson masses as a function of temperature. Also shown are 2m q and m q + m s . 
The Mott transitions of the pion, the kaon and the eta are marked by the arrows. 

FIG. 3. Effective quark masses (left) and temperatures (right) as a function of r for 
t = 0, 2 and 4fm/c. In the left panel, the upper curves are for strange quarks and the 
lower ones for light quarks, respectively. 

FIG. 4. Time dependence of the particle multiplicities. 

FIG. 5. Angular averaged particle density at t = lOfm/c. 

FIG. 6. Histogram of the creation temperatures of pions (solid) line, kaons (dashed 
line) and etas (dot-dashed line). 

FIG. 7. Particle density along the coordinate axes. Upper panel: t = 0, lower panel: 
t = 30fm/c. 
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